Parametric amplification of matter waves in dipolar spinor Bose-Einstein condensates 
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Spin-changing collisions may lead under proper conditions to the parametric amplification of matter waves in 
spinor Bose-Einstein condensates. Magnetic dipole-dipole interactions, although typically very weak in alkaline 
atoms, are shown to play a very relevant role in the amplification process. We show that these interactions may 
lead to a strong dependence of the amplification dynamics on the angle between the trap axis and the magnetic- 
field orientation. We analyze as well the important role played by magnetic-field gradients, which modify 
also strongly the amplification process. Magnetic-field gradients must be hence carefully controlled in future 
experiments, in order to observe clearly the effects of the dipolar interactions in the amplification dynamics. 
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I. INTRODUCTION 

Spinor Bose-Einstein condensates (BECs), formed by 
atoms with various available Zeeman states, have attracted a 
large attention in recent years, mostly motivated by the rich 
physics resulting from the interplay between internal and ex- 
ternal degrees of freedom. In addition to a wealth of possi- 
ble ground-state phases [1-4], the spinor dynamics has been 
at the focus of major interest ||5l |a]. This dynamics results 
from spin-changing collisions, which coherently re-distribute 
the populations among the different Zeeman states. Interest- 
ingly, spin-changing collisions are typically characterized by a 
very low energy scale much lower than the chemical potential 
in the condensate. As a result of that, the spinor dynamics in 
alkaline gases may be extraordinarily sensitive to other small 
energy scales. 

Up to very recently, only short-range interactions have 
played a role in typical experiments in ultra-cold gases. Re- 
cent experiments have started to unveil the rich physics re- 
sulting from the dipole-dipole interactions (DDI) IJSl |9|] . This 
is particularly the case of Chromium, which presents a rel- 
atively large magnetic dipole moment, /i = 6/iB. Remark- 
able effects of the magnetic DDI have been reported in recent 
experiments on Chromium BEC |10-13]. Alkaline atoms, 
on the contrary, present a much lower magnetic dipole mo- 
ment, /i = /^s/2, and hence they are not usually expected 
to show any trace of the DDI unless short-range interactions 
are switched-off by means of Feshbach resonances [ 14l[l5ll . 
However, as mentioned above, the spin-changing collisions in 
alkaline spinor BECs (in particular F = 1 *^Rb) are remark- 
ably low-energetic. As a result, spinor dynamics is very sensi- 
tive to magnetic DDI, in spite of the very low magnetic dipole 
moment. Recent experiments |16] have shown that the DDI 
may induce magnetization patterns in F = 1 *^Rb BECs. 

Spinor dynamics is particularly interesting for the case of 
condensates initially prepared in the tti = Zeeman sublevel. 
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In that case, spin-changing collisions may lead to correlated 
Einstein-Podolsky-Rosen (EPR) pairs in to = ±1 lIlTl [iST . 
in a process which closely resembles parametric down con- 
version in non-linear optics il9ll . As a result, spinor conden- 
sates may act as parametric amplifiers of matter waves 11201 - 
I22I1 . opening interesting perspectives for the creation of non- 
classical states of matter based on spinor BECs. Recently 
we have shown that the interplay between trapping potential, 
quadratic Zeeman effect (QZE) and spin-changing collisions 
crucially determines the amplification gain [21]] and its sensi- 
tivity with respect to quantum spin fluctuations 112211 . 

In this paper we show that the amplification dynamics may 
be extremely sensitive to the DDI. As a result of that, the am- 
plification of EPR-like pairs is strongly modified by the rel- 
ative orientation between the applied magnetic field and the 
trap axis. We analyze in detail this dependence, as well as the 
effects of magnetic-field gradients. We show that these gradi- 
ents must be carefully controlled, since uncontrolled gradients 
may obscure the expected DDI effects. We finally comment 
on experimental requirements. 

The structure of the paper is as follows. In Sec. |II| we 
present the system considered and the corresponding Hamil- 
tonian. The linear regime is discussed in Sec.|lll| An intuitive 
qualitative picture of the effects of the DDI in the amplifica- 
tion dynamics is discussed in Sec.|IV] In Sec. [V] we introduce 
the main formalism to analyze the amplification dynamics in 
the presence of DDI, whereas the coiTesponding numerical re- 
sults are presented in Sec. [Vlj The effects of the magnetic- 
field gradient are analyzed in Sec. IVIII Finally we discuss 
experimental requirements and summarize our conclusions in 
Ssc.NM 



II. HAMILTONIAN 



In the following we consider a spin-1 Bose gas (e.g. F = 1 
*^Rb), with Zeeman components m = 0,±1, confined in 
a dipole trap in the presence of an external magnetic field 
(which we assume as oriented along the z-axis). The system 
is described by the Hamiltonian H = Hq + H^r + ^dd- In this 



Hamiltonian, the single-particle physics is described by 



H^^j d\Y.^l{r) 



2M 



+ V{r) + Ez{m) 



i>mir), (1) 



where V'm annihilates bosons with spin projection m. The 
trapping potential is of the form 1^ (r) = 



{u^lia 



y 



\ 



where M is the atomic mass, cjii ^ lu± are the trap 



frequencies (cigar-shaped trap), and y' — cos dy + sin -d z, 
z' = — sin "dy + cos -d z, with d the angle between the trap 
axis and the magnetic field orientation. This angle will play a 
crucial role in our discussion of the effects of the DDL 
The Zeeman energy for the m component is of the form: 



Ez{m) ~ (p + Vp ■ f) rn + qm , 



(2) 



where p — gLfJ-BBo characterizes the linear Zeeman ef- 
fect (LZE) for a homogeneous magnetic field Bq, with g^ 
the Lande factor {gL = -1/2 for F = 1 ^'^Rb) and ^s 
the Bohr magneton. The quadratic Zeeman effect (QZE) is 
characterized by a constant q, which in principle depends as 
q = fi^ Bq / (SChfs) on the hyperfine coupling strength Chts 
(w ft, X 3.4 GHz for ^^Rb), but may be also externally modi- 
fied using optical or micro-wave dressing ll23ll24|] . Addition- 
ally, we allow for a magnetic field gradient B = Bo + VB ■ r, 
leading to an energy shift m Vp • r, which plays a relevant role 
below. 

The short-range interactions are given by: 



H, 



d'rY, V'i.(r)Vi,L.(r)C/r.r^^'(r)^«(rO, (3) 



with U'^;^ = UoSm..fhSm',rh' + Uifmrh ' fm'rn', wherC 

/mm' = ifLn'Jlm'JLn'V^ with f^^y^, the spin-1 ma- 
trices. C/q ~ (go + 2.g2)/3 and Ui = {g2 — .9o)/3 are, re- 
spectively, the coupling constants for spin-independent and 
spin-dependent interactions, where gp — Anh'^aF/M, with 
ap the s-wave scattering length for the channel with total spin 
F. Note that the short-range interactions preserve the total 
spin projection, but this may be done in two crucially dif- 
ferent ways, either by preserving the individual spin projec- 
tions (spin-preserving collisions) or by modifying the individ- 
ual projections while preserving the total one (spin-changing 
collisions). 

Finally, the magnetic dipole-dipole interaction is given by: 
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dd 
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with 



W , (r — r I = 
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where (P — /io(7£/i^/(47r), /Uq is the vacuum permeability, 
and Ur — [r ~ r')/\r — r'\. Contrary to the short-range in- 
teractions the DDI may violate, in principle, the conservation 
of the total spin projection ( they may induce the equivalent of 
the Einstein-de Haas effect ll25ll26ll ). However, the associated 
change in LZE is typically, even for very low magnetic fields, 
orders of magnitude larger than any energy in the system and 
hence these spin-violating processes can be safely considered 
as suppressed. 

Since short-range interactions preserve spin projection, and 
so do in practice DDI as well, the homogeneous LZE {pm) 
is preserved and it may be gauged out, playing no role in the 
dynamics discussed below. The same argument cannot be, 
however, employed with the magnetic-field gradient which 
may play a significant role in the spinor dynamics 1127(1 . as 
discussed in Sec.lVllI 



III. LINEAR REGIME 

In the following we are interested in the first stages (lin- 
ear regime) of the spinor dynamics of a spin-1 BEC initially 
prepared in the to = sublevel, after quenching q into the un- 
stable regime. This dynamics, induced by spin-changing col- 
lisions, is characterized by the correlated creation of atomic 
pairs in TO, = ±1. In this section, and for the sake of simplic- 
ity, we do not consider magnetic-field gradients, which will 
be introduced in Sec. lVIII 

Before quenching q, the BEC is prepared in the to = 
component. The initial scalar wavefunction i/;o of the BEC in 
TO = and the corresponding chemical potential jjl may be ob- 
tained from the time-independent Gross-Pitaevskii equation: 



2M 



A + y(r) + C/ono(f) 



V'oC'^) =M'0o(?^), (6) 



withno(f) = |^o(^)p, snA j (Pr\%l)Q{r)\'^ = TV. 

The first stages of the spinor dynamics may be described by 
means of a Bogoliubov approximation: 



(Vi, V'o, V'-i)^ = (0, ^0, 0)^ + Wi, <5^o, &i>~i? 



-ifit 



(7) 



where we consider small fluctuations of the spinor field oper- 
ator {Slpni}, such that |V'Op > Y.mi^'^LS'^m)- 



A. Hamiltonian without dipole-dipole interactions 

We consider first the Hamiltonian without DDL Insert- 
ing (|7]i into the grand canonical potential Hq + Hsr — fJ-N 
(with N — J (PrJ2m.''Pmi''^)^rn{r) the total particle num- 
ber), and keeping terms up to second order in Stpm, we obtain 
an effective Hamiltonian for Sip±i of the form: 



^1 = E fd'rS^l{H,ff + q)Si,ra 

m=±l-' 

+ Ui [ d\ no {6^pl6'il>U + HM-i) ■ (8) 



Note that in the linear regime the fluctuations Sip±i are de- 
coupled from the density and phase fluctuations 5ipo of the 
m = BEC (which may be excited during the prepara- 
tion process). In the previous expression we have introduced 

H^ff = -h'^A/2M + Veff{f), where 



^eff 



K//(r) - V{r) + {Uo + Ui)noir} - fi 



(9) 



may be understood as the effective trapping potential felt by 
the ±1 -fluctuations on top of the to = BEC. It contains the 
mean-field potential {Uo + Ui)na, which originates from spin- 
preserving collisions of ±l-atoms with the BEC in m = 0. 
Note that in the Thomas-Fermi regime fi — ¥{7^) + Uono{f). 
In that regime, Ve//(r) = Uino{f) within the BEC region, 
and Veff = V{r) — fi outside. 

The second line in Eq. ^ originates from spin-changing 
collisions, which convert to = atoms into ±l-atom pairs 
and vice-versa. Interestingly, this process resembles paramet- 
ric down conversion in optical parametric amplifiers [|l9tl . In- 
deed, if the m — BEC is unstable after the quench of q, 
spin-changing collisions lead to an exponential amplification 
of the population in to = ±1 0201 - 12211 . 



B. Dipole-dipole interactions 

As mentioned above, spin-changing collisions are typically 
characterized by a very low energy scale. As a result, the 
spinor physics is highly sensitive to other very low energy 
scales, including the rather weak magnetic DDI in alkaline 
gases. This sensitivity has been recently demonstrated in 
experiments on the formation of spatial magnetization pat- 
terns in spinor Rb BECs [16]. As we show in following sec- 
tions, also the exponential amplification of the population in 
771 = ±1 following a quench in q may be very significantly 
modified by the DDI. 

Although there is no DDI contribution in the GP Eq. (|6|, 
since /oo = 0, there is however an important contribution 
to the effective Hamiltonian for 6t/j±i which we obtain after 
inserting (|7]l in Hdd and linearizing: 



Hi 



dd 



d'r j d\'Mr)Mr')Vii{r-f') 



+ 5i^\{f)5i^l^(r')+5il:i{r)5i}-i{f') 



(10) 



with Vdd(r) == ^ (3z2 - |f|2). In Eq. ^ we have ne- 
glected terms related to scattering processes which do not pre- 
serve the total spin projection since, as mentioned in Sec. [Ill 
the associated change in LZE suppresses spin-violating pro- 
cesses even for very low magnetic fields. Note that the third 
line of Hidd contains as well a parametric amplification term. 



IV. QUALITATIVE PICTURE OF THE EFFECT OF THE 

DIPOLE-DIPOLE INTERACTIONS ON THE 

AMPLIFICATION DYNAMICS 

The effects of the DDI on the amplification dynamics may 
be qualitatively understood from a simplified model. In homo- 
geneous space (V — 0, constant tiq, /^ = UqUq) we may in- 
troduce the Fourier transform fim{k) ~ J cPr 5'iprn{'r)e~'^^'^ , 
which allows to write Hi in the simpUfied form ll28ll : 
i/^^o™ = / hkd^k/{2TT)'\ where 

hk = Yl i^k + q - qcr) vLWVrnik) 



m=±l 



- qcr [4iM-ii-k) + mik)fl-i{-k)j , (11) 

with Ek ^ h?k'^ /2M and q^r = -UiUq. Note that for F = 1 
*^Rb Ui < 0, and hence qcr > (we consider in the following 
this case, although for the F — 2 case the sign is the opposite). 
This Hamiltonian, which may be easily diagonalized for each 
momentum fc, possesses eigenenergies of the form 



\^{k) = ViEk+q-qcrY 



(12) 



Note that if Im(A^ (fc)) > for some eigenenergy, then there 
is an exponential growth of spin excitations, which leads to a 
correlated pair creation in to = ±1. This instability is best 
characterized by the instability rate A ~ max{Im(A±(fc))}. 
It is straightforward to show that the to = BEC becomes 
unstable for q < 2qcr- The instability rate rises between 
qcr < q < 2qcr acquiring its maximal value A = q^r at 
q — qcr- For q < qcr, A = qcr in the homogeneous case 
(in the inhomogeneously trapped case the instability rate A 
presents significant modulations which are responsible of the 
multi-resonant g-dependence of the amplification dynamics 
recently observed experimentally 1I21II ). 

In the following we apply a similar formalism to the DDI 
term Hidd- We introduce the Fourier transformation of 
Vdd{r) to obtain Vdd{k) = C/drf(l - 3cos^ 9), where Udd = 
2TTd? /?> and is the angle between k and the dipole orientation 
(z axis). Using convolution theorem, and since we assume tiq 
as constant, we may then rewrite Hi^dd in the form: 



Tjhom 
^l,dd 



(Fk 

W) 



7 no 



Vddik)\ J2 fll{k)flm{k) 



,m=±l 



r}l(fc),7_i(-fc)t+,}i(fc),}_i(-fc)} 



(13) 



Note that adding Hi^cid to Hi results in a similar form as that 
of Hi but with an effective q^l^ {9) — qcr — noVdrf(fc). Note 
that, due to the anisotropy of the DDI, q'^cV depends on the 
angle 9. The effects of the trap geometry may be qualitatively 
understood from this 9 dependence. For an axisymmetric trap 
the dominant momenta are those along the tightest direction. 
If the dipole orientation is perpendicular to the trap axis, then 
the dominant k will then be those with = 0, and hence 
q^cP — Qcr + '^iT'oUdd- On the contrary, if the dipole orien- 
tation is parallel to the trap axis, the dominant momenta will 
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FIG. 1 : (Color online) Instability rate A as a function of g for A'^ = 
10^CJ± = 27rx200Hzanda;|| = 27r x 40 Hz, for the case of no DDI 
(black, dashed), Q — n/2 (blue, solid) and ^ = (red, dot-dashed). 



be those with 6 = 7r/2, and q[ 



eff 



qc 



noUdd- Since 



A ~ Qcr ' ^^ hence expect an enhancement of the instabil- 
ity for a magnetic field orientation perpendicular to the trap 
axis, and a reduced instability for a parallel orientation. Al- 
though the DDI in alkaline atoms is typically very weak, the 
spin-changing collisions are very weak as well. In particu- 
lar, in i^ — 1 ^^Rb the strength of the DDI is quite signifi- 
cant compared to the strength of the spin-changing collisions, 
\Udd/Ui\ w 0.2. As a result, the DDI modification of the 
instability rate is expected to lead to a marked orientation de- 
pendence of the amplification dynamics. In the following sec- 
tions we show that this is indeed the case when considering 
realistic trapped cases. 



V. AMPLIFICATION DYNAMICS IN TRAPPED DIPOLAR 
CONDENSATES 

Although the homogeneous model discussed before allows 
for a simplified intuitive understanding of the major effects 
of the DDI in the amplification process, a quantitative anal- 
ysis of realistic experimental situations may be just achieved 
by properly considering the inhomogeneous trapping, and the 
corresponding inhomogeneous density no (r) of the m = 
BEC. In this section we introduce the basic formalism which 
we follow for the analysis of the amplification dynamics char- 
acterizing the spinor physics in the linear regime. 

The analysis of the spinor dynamics is significantly sim- 
plified by considering the eigenfunctions and eigenener- 
gies of Heff, H^ff(f>n{f) = e„0„(r), and expanding the 
field operators in the basis of these eigenstates 6iprn,{r) ~ 
J2n 0n(^)«nm- We may then rewrite: 



Hi+Hidd = 



/ ^vv^n ~r QjOfiji' -r iDjui' ) ^ ^ ^nrn^Ti'm 



^t 



m=±l 



+ y^^{Ann' + Bnn') (alia|,/_i + anlSn'-lj , (14) 



where A„„/ = Ui J d^r no4>n4'n' characterizes the effects of 
the short-range spin-changing collisions, whereas the effects 
of the DDI are given by 

Bnn'^fd^ J d''r'Fr,{f)Vddir- f )^«'(^), (15) 

where Fn{r) — '0o(r)0n(?O- The matrix elements Bnn' are 
most efficiently calculated in fc-space according to 

d'^k ^ ^ ^ 

Bnr,- = I J^Fnik)Vddik)Fn'{k) 



(16) 



where F{k) is the Fourier transform of F{'r). 

Eq. ( fT4] i is solved by the multimode Bogoliubov ansatz 



E 






where d^ satisfy a J , i/i 
to the eigenvalue equation: 



Hi. 



dd 



Xtat, 



(17) 



which leads 



C 



-± 



A 



(18) 



where u 



■±T — 



^b 



E + 91 + 

A + B 



, } (and similarly for v^) and 

B ~A + B) 
-E 9I-B 



(19) 



with Enn' — ^nSnn', 1 the identity matrix, and A (B) the 
matrix with components Ann' {Bnn')- From the Heisenberg 
equations of motion a^{t) — Q;J(0)e~*^'' */^. Note that, as 
for the homogeneous case, if Im(A^) > for some eigenen- 
ergy, then there is an exponential growth of correlated pairs in 
771 = ±1. As for the homogeneous case, this instability is best 
characterized by the instability rate A ~ max{Im(Ai^)}. 
The time evolution of a„_±i is then easily obtained 



{ani{t)} 



{aLiWI 



U(t) 



{a„i(0)}' 
{SLi(O)}. 



(20) 



with U — M^^e~*'^*/''M, where M is the matrix of eigen- 
vectors obtained after solving Eq. ( fTSl l and A the correspond- 
ing diagonal matrix of eigenvalues. 

As mentioned above, the atoms are initially prepared in the 
7n = sublevel. However, a slightly imperfect preparation 
may lead to a non-zero population of Ng m = ±.1 pairs in the 
original BEC [22.1 . These spurious atoms (which from now are 
called classical seed) share the same wavefunction ijjoi'r) as 
the 771 = atoms. Denoting Xn = f d'^ripa{r)<j)nir)/VN, we 



^n,m'-^n,ml 



may then easily express the population P^ — J2r 
in the form P„(t) = Pc{t) + Pgit), where 

Pc{t) = Nsx-(o^o + 6^6)x, (21) 

denotes the population triggered by the classical seed, and 



PQ{t)= Tr[0^0 



^t, 



(22) 



denotes the population induced by quantum fluctuations (i.e. 
when Ng = 0). In the previous expressions, the matrices O 
and O are the upper left and upper right part of the time evo- 
lution matrix U(t) and X = {xi,X2, ■■■)'^ ■ 




4 6 

g//*(Hz) 



FIG. 2: (Color online) Fraction of atoms transfered into | ± 1) after 
115 ms as a function of q for Na — 2 and the same parameters as 
Fig.m and for the case of no DDI (black, dashed), 6 — -k 12 (blue, 
solid) and 6 = (red, dot-dashed). 
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FIG. 3: (Color online) Maximal population transfer after 115 ms as 
a function of Q (same parameters as Fig. [2} for different magnetic 
field gradients VB = (blue, solid), 5 mG/cm (black, dashed) and 
lOmG/cm (red, dot-dashed). 



VI. DIPOLE-INDUCED ORIENTATION-DEPENDENCE OF 
THE AMPLIFICATION DYNAMICS 

In this section we employ the formalism discussed in Sec.lVl 
to study the effects of the DDI in the amplification dynamics. 
We shall show that due to the DDI the amplification may be 
markedly dependent on the relative orientation between the 
trap axis and the external magnetic field. 

In our numerical calculations we have considered realistic 
experimental conditions, with iV = 10^ i^ = 1 ^^Rb atoms in 
a cigar-shaped harmonic potential with wj^ = 27r x 200 Hz, 
and wy = 27r X 40 Hz. As mentioned above, we consider 
the atoms as initially prepared in m = with possibly an 
initial spurious classical seed (which we typically consider as 
A^s = 2, a typical value expected from previous experimental 
results 1221 ). At i = the QZE energy q is set to a given value 
within the instabiUty regime. We monitor the subsequent evo- 
lution of the populations P±i(t) obtained from Eqs. (|2TI) and 
(|22] | as a function of q and the relative angle 6 between the 
trap axis and the external magnetic field. 

Fig. [T] shows the dependence of the instability rate A as a 
function of q for 6* = 0, 6* = 7r/2 and without DDI. Note that 
for all cases, the instability rate experiences a maximum con- 
trary to the homogeneous case. This maximum is induced by 
the inhomogeneous harmonic trapping and leads to marked 
resonances in the g-dependence of the amplification dynam- 
ics, as discussed in Ref. 1.2 1.1 . However, A clearly depends on 
the trap orientation confirming indeed the intuitive qualitative 
picture discussed in Sec. HV] When trap axis and magnetic 
field are aligned A decreases compared to the non-dipolar 
case, whereas the opposite is true when the magnetic field is 
oriented perpendicular to the trap axis. Note as well that, also 
as expected from the qualitative picture of Sec. IIVI the insta- 
bility region is shifted towards lower q values in the parallel 
configuration, and towards larger q values in the perpendicular 
one. 

This modified instability rate translates into a significantly 
distorted pair-creation dynamics, due to the exponential na- 



ture of the parametric amplification. Fig. |2] shows the trans- 
ferred fraction P±i{t)/N after t = 115 ms as a function of 
q for different values of 9. As expected from the form of A 
we observe the appearance of a maximum for all 6, which 
is slightly shifted (by approximately 1 Hz) towards lower q 
when 9 is shifted from 7r/2 to 0. However, this maximum is 
approximately four times as large for 9 = 7r/2 than for 6^0. 
The dependence of the amplification on 9 is very clearly ob- 
servable in the 6'-dependence of the maximum of P±i (again 
at i = 115 ms) shown in Fig.|3] Note that the maximum grows 
monotonically from 6* = to 6* = n/2. 



VII. EFFECTS OF MAGNETIC FIELD GRADIENTS ON 
THE AMPLIFICATION DYNAMICS 

As mentioned in Sec. |II]the homogeneous LZE plays typ- 
ically no role in the spinor dynamics (only at very low mag- 
netic fields B < 1 mG the DDI could induce the equivalent 
of the Einstein-de Haas effect ll25ll26ll . and in this case the 
residual LZE could play a role). However, magnetic field gra- 
dients cannot be gauged out, and may play a relevant role in 
the spinor physics 1271 . In this section, we analyze the effects 
that these gradients may have on the amplification dynamics. 
We shall show that even relatively weak gradients may have a 
significant effect on the amplification process. 

Although magnetic-field gradients do not affect the GP 
equation for the rn — BEC, there is indeed a contribution 
to the effective Hamiltonian for (5V'±i in the linear regime: 



Hl,gr = Vp • 



d-^r 



(si^lrS^jji - Sipl^fSijj-A , (23) 



which may be straightforwardly implemented into the matrix 
C of the eigenvalue equation ( fTSl l, 



C = 

with D„ 



E + ql + B + B -A + B) 

A + B Egl-B + D 

= Vp- J d^r(f>nr(l)n'. 



(24) 
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FIG. 4: (Color online) Fraction of atoms in ±1 after a time 115 ms 
(same parameters as Fig. |2]l as a function of g for ^ — tt/2, and 
for gradients VB — (blue, solid), 5 mG/cm (black, dashed) and 
10 mG/cm (red, dot-dashed). 



The magnetic-field gradients have two main effects. On one 
side, they modify the effective potential Vg// (r) in a different 
way for m — I than for m ~ —1. This reduces the overlap of 
the 771 = ±1 atom clouds with the m = BEC and hence 
the scattering mediated transfer On the other side, atoms 
placed at different locations experience different Larmor pre- 
cession frequencies. Although this does not affect the local 
short-range interactions, it does modify the non-local DDL 
For large-enough gradients this may lead to a time-averaged 
DDI 1 27]. For weak gradients, as those considered below, the 
explicit time dependence induced by the gradients must be 
considered. 

Parametric amplification is handicapped by the presence of 
gradients as a result of these two combined effects. Fig. |4] 
shows the combined effect of the DDI and the magnetic-field 
gradient along the weak trap axis for 9 ~ 7r/2. As expected, 
we obtain a reduction of the transfer maximum with increas- 
ing gradient and a shift of its position to lower q/h by approxi- 
mately 1.5 Hz for a gradient of 10 mG/cm. Hence the transfer 
maximum is shifted down and to lower q with decreasing 9 
and increasing gradient. As shown in Fig. [3] in the presence 
of a magnetic field gradient the maximum of P±i{t) shows 
also a marked 9 dependence. 

Hence, even rather weak gradients (< 10 mG/cm) may 
strongly modify the amplification dynamics, an effect which is 
enhanced by the presence of the DDI. Although as mentioned 
above, the 0-dependence should reveal also in the presence 
of gradients the effects of the DDI, slight variations of the 
magnetic field gradients (of the order of few mG/cm) when 
changing the magnetic-field orientation with respect to the 
trap axis must be very carefully controlled. This is indeed 
a crucial point, since otherwise, accidental ^-dependences 
of the magnetic-field gradients may obscure the physical 9- 
dependence characteristic of the DDI. 



VIII. DISCUSSION 
A. Experimental requirements 

In the following we outline the general requirements to per- 
form an experiment with F = 1 ^^Rb to probe the theoretical 
results discussed above. A m = BEC must be prepared in a 
crossed dipole trap following the same procedure detailed in 
Refs. ||2U|22I1 . The state preparation requires particular atten- 
tion, since remaining atoms in m 7^ may strongly alter the 
experimental result. Previous experiments have shown that 
the number of atoms in m 7^ states can be suppressed to 
TVs w 2 by briefly applying a strong magnetic field gradient to 
purify the system [22]. Due to the nature of the QZE in F = 1 
^^Rb, magnetic fields can be used to access positive values of 
q, as required above. 

As discussed in Sec. IVIII the most significant requirement 
compared to previous experiments is related to the suppres- 
sion of magnetic field gradients which could obscure the dipo- 
lar effects. Appropriate experiments should be carefully de- 
signed to minimize all sources of field gradients from the 
vicinity of the atomic sample (alternatively a magnetic shield 
could be placed around the sample). In state of the art preci- 
sion measurements, field gradients are commonly suppressed 
below 1 mG/cm [29], which is sufficient to realize the manda- 
tory experimental conditions (see Sec. lVIII l. 

Fig.|2]shows that the resonant spin transfer to the m = ±1 
states depends strongly on the relative orientation of the weak 
trap axis and the external magnetic field. Since it is difficult 
to change the orientation of a dipole trap while maintaining 
its trapping potential, experiments must be designed to vary 
the orientation of the external magnetic field. In this sense, 
two sets of Helmholtz coils are necessary to provide a ho- 
mogeneous external magnetic field. One of them should be 
placed along the weak axis of the trapping potential to real- 
ize the 9 — configuration and another one along one of 
the strong axis to realize the 9 — tt/2 case. Both mag- 
netic fields have to be calibrated, preferentially using preci- 
sion microwave spectroscopy between the ground state hy- 
perfine manifolds of *^Rb. Such an experimental apparatus 
would also allow for a rotation of the field, since the currents 
in the two Helmholtz coils could be adjusted to obtain a rel- 
ative angle 6. In this way, it should be possible to perform a 
measurement analogous to that discussed in Fig. [3] Finally, 
additional magnetic field gradients can be applied along both 
magnetic field directions to observe the suppression shown in 

Fig. a 



B. Summary 

We have shown that, in spite of the very small magnetic 
moment, the magnetic DDI may lead to a strong modifica- 
tion of the amplification dynamics in F = 1 *^Rb due to the 
low-energy scale of the spin-changing collisions. We have 
shown that the DDI induce a very marked dependence of the 
amplification gain with respect to the relative orientation be- 
tween magnetic-field direction and trap axis. If both directions 



are perpendicular to each other the amphfication dynamics is 
much faster than for the parallel configuration. Remarkably, 
the number of transferred atoms into m ^ ±1 may increase 
for F — I ^^Rb for a fixed holding time of around 100 ms by 
a factor over 400% when turning from a parallel to a perpen- 
dicular configuration. We have shown as well that magnetic 
field gradients may also significantly modify the amplification 
dynamics, both due to their effects on the trapping and on the 
DDL As a result, magnetic-field gradients must be carefully 
controlled, since uncontrolled changes in the gradient when 
turning the magnetic field orientation may obscure the orienta- 



tion dependence of the DDI effects on the amplification. This 
demands specific requirements for future experiments. 
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